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I. INTRODUCTION 

The production of gluons and quarks in high-energy 
coUisions, and their development into the jets of parti- 
cles observed in experiments, is usually described by the 
theory of Quantum Chromodynamics (QCD). In pertur- 
bative QCD, a produced parton (gluon or quark) emits 
gluon radiation, with each subsequent emission carrying 
off a fraction of the original parton's energy and mo- 
mentum. The probability for a gluon to radiate a gluon 
is proportional to the color factor Ca = 3, while gluon 
radiation from a quark is proportional to the color fac- 
tor Cp = 4/3. In the asymptotic limit, in which the 
radiated gluons carry a small fraction of the original par- 
ton's momentum, and neglecting the splitting of gluons 
to quark-antiquark pairs (whose probability is propor- 
tional to the color factor Tr — 1/2), the average number 
of objects radiated by a gluon is expected to be a factor 
Ca/Cf — 9/4 higher than the number of objects radi- 
ated by a quark 0. In general, it is expected that a 
gluon will yield more particles with a softer momentum 
distribution, relative to a quark |^,^ . 

Although gluon jets are expected to dominate the final 
state of proton-antiproton (pp) collisions at high energies, 
quark jets make up a significant fraction of the jet cross 
section at high xt — ^=-, where ^/s is the total energy 

of the pp system, and pt is the jet momentum transverse 
to the hadron-beam direction. The ability to distinguish 
gluon jets from quark jets would provide a powerful tool 
in the study of hadron-collider physics. To date, however, 
there has been only little experimental verification that 
gluon jets produced in hadron collisions display charac- 
teristics different from quark jets For fixed pr, 
we analyze the internal structure of jets at y/s = 1800 
GeV and 630 GeV by resolving jets within jets {sub- 
jets) RPHlTt- Using the expected fractions of gluon and 
quark jets at each y/s, we measure the multiplicity of sub- 
jets in gluon and in quark jets. The results are presented 
as a ratio of average multiplicities r = of subjets 
in gluon jets to quark jets. This measured ratio is com- 
pared to that observed in e+e~ annihilations J^,^] , to 
predictions of a resummed calculation ||ll| , [l^ , p7[j 7 and to 
the HERWIG jl^ Monte Carlo generator of jet events. 

The D0 detector |l|], described briefiy in Sec. ||, is 
well-suited to studying properties of jets. A jet algo- 
rithm associates the large number of particles produced 
in a hard-scattering process with the quarks and gluons of 
QCD. We define jets with a successive-combination algo- 
rithm p0|-p3t based on relative transverse momenta (fcj^) 
of energy clusters, described in Sec. III. In this paper, 
we present the first measurement of jet properties using 
the k± (sometimes written kx) algorithm at a hadron 
collider. The momentum calibration of jets in the fcj^ 
algorithm is outlined in Sec. Ill C| , followed in Sec 



scribed in Sec. [HE. In e+e annihilations, the number 



of subjets in gluon jets was shown to be larger than in 
quark jets |13 1^. In pp collisions, identifying gluon and 
quark jets is more complicated than in e~^e~ annihila- 
tions. We approach this issue by comparing central jet 
samples at y/s — 1800 GeV and 630 GeV, with the sam- 
ples described in Sec. IV. For moderate jet pr (55 - 
100 GeV), the y/s = 1800 GeV sample is gluon-enriched, 
and the y /s = 630 GeV sample is quark-enriched. Sec- 
tion IV D| describes a simple method developed to extract 
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by a simple comparison with jets defined with the fixed- 
cone algorithm. To study jet structure, the k± algorithm 
is then applied within the jet to resolve subjets, as de- 



thc separate subjet multiplicity for gluon and for quark 
jets. The method does not tag individual jets, but in- 
stead, we perform a statistical analysis of the samples at 
y/^ = 1800 GeV and 630 GeV ||. The method requires 
the relative mix of quarks and gluons in the two data 
samples, which is derived from a Monte Carlo eve nt gen - 
erator that uses the parton distribution functions ||25| , p6| , 
measured primarily in deep inelastic scattering. Subse- 
quent sections describe the measurement of the subjet 
multiplicity in D0 data and Monte Carlo simulations, 
the corrections used in the procedure, and the sources of 
systematic uncertainty. We conclude with comparisons 
to previous experimental and theoretical studies. 



II. D0 DETECTOR 

D0 is a multipurpose detector designed to study pp 
collisions at the Fermilab Tevatron Collider. A full de- 
scription of the D0 detector can be found in Ref. [p^ . 
The primary detector components for jet measurements 
at D0 are the excellent compensating calorimeters. The 
D0 calorimeters use liquid-argon as the active medium to 
sample the ionization energy produced in electromagnetic 
and hadronic showers. The elements of the calorime- 
ter systems are housed in three cryostats. The cen- 
tral calorimeter (CC) covers the region I77I < 1.0, while 
the symmetric end calorimeters (EC) extend coverage to 
|ry| < 4.2, where the pseudorapidity 77 — — lntan6'/2 is 
defined in terms of the polar angle 9 with respect to the 
proton-beam direction z. Each system is divided into 
an electromagnetic (EM), fine hadronic (FH), and coarse 
hadronic (CH) sections. The EM and FH use uranium 
absorber plates as the passive medium, and the CH uses 
either copper (CC) or stainless steel (EC). Copper read- 
out pads are centered in the liquid-argon gaps between 
the absorber plates. Radially, the electromagnetic sec- 
tions are 21 radiation lengths deep, divided into 4 read- 
out layers. The hadronic calorimeters are 7-11 nuclear 
interaction lengths deep, with up to 4 layers. The en- 
tire calorimeter is segmented into towers, of typical size 
A77 X Acf) = 0.1 X 0.1, projected towards the nominal pp 
interaction point in the center of the detector, where 
is the azimuthal angle about the z axis. Figure |l| shows 
a schematic view of one quadrant of the D0 calorime- 
ter in the r — z plane, where r is the distance from the 
origin in the plane transverse to the beam axis. Each 
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FIG. 1. One quadrant of the D0 calorimeter and drift 
cliambers, projected in the r — z plane. Radial lines illus- 
trate the detector pseudorapidity and the pseudo-projective 
geometry of the calorimeter towers. Each tower has size 
A»7 X A<^ = 0.1 X 0.1. 



layer in a calorimeter tower is called a cell, and yields 
an individual energy sampling. Energy deposited in the 
calorimeters by particles from pp collisions are used to 
reconstruct jets. The transverse energy resolution of jets 
for data at y/s — 1800 GeV can be parameterized as [ETl: 



{a{ET)/ET)' 



+ 0.5 /Et + 0.001, (2.1) 



with Et in GeV. 

In the analysis of jet structure, we are interested in 
the distribution of energy within jets. Apart from the 
energy of particles produced in a hard-scattering event, 
the cells of the D0 calorimeter are sensitive to three ad- 
ditional sources of energy that contribute to a jet. The 
first, called uranium noise, is a property of the detec- 
tor material. The decay of radioactive uranium nuclei 
in the calorimeter can produce energy in a given cell, 
even in the absence of a particle flux. For each cell, 
a distribution of this pedestal energy is measured in a 
series of calibration runs without beams in the acceler- 
ator. The pedestal distribution due to uranium noise is 
asymmetric, with a longer high-end tail, as illustrated in 
Fig. 1^. During normal data-taking, the mean pedestal 
energy is subtracted online from the energy measured in 
a hard-scattering event. To save processing time and re- 
duce the event size, a zero-suppression circuit is used, 
whereby cells containing energy within a symmetric win- 
dow about the mean pedestal count are not read out. 
Since the pedestal distribution of each cell is asymmetric, 
zero-suppression causes upward fluctuations in measured 
cell energies more often than downward fluctuations. In 
the measurement of a hard-scattering event, the net im- 
pact is an increased multiplicity of readout cells and a 
positive offset to their initial energies. 

There are two other environmental effects that con- 
tribute to the energy offset of calorimeter cells. The 
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FIG. 2. Illustration of the pedestal energy distribution in 
a calorimeter cell (solid line), stemming from uranium noise. 
The mean value is defined to be zero, and the peak occurs at 
negative values. Removal of the portion between the vertical 
dashed lines (a symmetric window about the mean) yields a 
positive mean for the remaining distribution. 



first is extra energy from multiple pp interactions in the 
same accelerator-bunch crossing, and this depends on 
the instantaneous luminosity. To clarify the second ef- 
fect, called pile-up, we turn to how calorimeter cells are 
sampled, as is illustrated in Fig. |^. The maximum drift 
time for ionization electrons produced in the liquid-argon 
to reach the copper readout pad of a calorimeter cell is 
about 450 ns. The collected electrons produce an elec- 
tronic signal that is sampled at the time of the bunch 
crossing (base), and again 2.2 /is later (peak). The dif- 
ference in voltage between the two samples (peak relative 
to base) defines the initial energy count in a given cell. 
Because the signal fall-time (~ 30 /is) is longer than the 
accelerator bunch spacing (3.5 ^s), the base and peak 
voltages are measured with respect to a reference level 
that depends on previous bunch crossings. The signal 
from the current bunch crossing is therefore piled on top 
of the decaying signal from previous crossings. When 
a previous bunch crossing leaves energy in a particular 
cell, that cell's energy count will therefore be reduced on 
average, after the baseline subtraction. 
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FIG. 3. Schematic of signal voltage in a calorimeter cell 
as a function of time. The solid line represents the contri- 
bution for a given event (the "current" pp bunch crossing). 
In the absence of previous bunch crossings, the cell is sam- 
pled correctly at tt, just before a crossing, to establish a base 
voltage, and at tp, to establish a peak voltage. The voltage 
difference AV^ = V{tp) — V{tb) is proportional to the initial 
energy deposited in the cell. The dashed lines show exam- 
ple contributions from a previous bunch crossing containing 
three different numbers of pp interactions. The observed sig- 
nal is the sum of the signals from the current and previous 
crossings. (The figure is not to scale.) 



III. k± JET ALGORITHM 



The energy-weighted centroid of a cone is defined by: 



V 



(3.1) 



where the sum is over all particles i in the cone. The 
centroids are used iteratively as centers for new cones in 
[r], (j)) space. A jet axis is defined when a cone's cen- 
troid and geometric center coincide. The fixed-cone jet 
algorithm allows cones to overlap, and any single particle 
can belong to two or more jets. A second parameter, and 
additional steps, are needed to determine if overlapping 
cones should be split or merged [Q. 

The fcj^ jet algorithm offers several advantages over 
the fixed-cone jet algorithms, which are widely used 
at hadron colliders. Theoretically, the algorithm is 
infrared-safe and collinear-safe to all orders of calcula- 
tion 1 20 3^. The same algorithm can be applied to par- 
tons generated from fixed-order or resummation calcula- 
tions in QCD, particles in a Monte Carlo event generator, 
or tracks or energy depositions in a detector. 

The k± jet algorithm is specified in Sec. [II A. In 
Sec. IIIB, we describe the preclustering algorithm, the 



goal of which is to reduce the detector-dependent aspects 
of jet clustering (e.g., energy thresholds or calorimeter 
segmentation). The mom entum cali bratio n of fc_L jets 

jets recon- 



IIID 



is presented in Sec. |III C| . In Sec 
structed using the k±_ algorithm are compared to je ts re- 
constructed with the fixed-cone algorithm. In Sec. [II E 
we indicate how subjets are defined in the fc_L algorithm. 



Jet algorithms assign particles produced in high-energy 
collisions to jets. The particles correspond to observed 
energy depositions in a calorimeter, or to final state par- 
ticles generated in a Monte Carlo event. Typically, such 
objects are first organized into preclusters (defined be- 
low), before being processed through the jet algorithms: 
The jet algorithms therefore do not depend on the nature 
of the particles. We discuss two jet algorithms in this pa- 
per: the fcj^ and cone jet algorithms, with emphasis on 
the former. 

In the k±^ jet algorithm, pairs of particles are merged 
successively into jets, in an order corresponding to in- 
creasing relative transverse momentum. The algorithm 
contains a single parameter D (often called R in some ref- 
erences), which controls the cessation of merging. Every 
particle in the event is assigned to a single k±_ jet. 

In contrast, the fixed-cone algorithm associates 
into a jet all particles with trajectories within an area 
A — ttTZ^ , where the parameter TZ is the radius of a cone 
in {rj, (f)) space. The D0 fixed-cone algorithm 1^,^ is 
an iterative algorithm, starting with cones centered on 
the most energetic particles in the event (called seeds). 



A. Jet clustering 

There are several variants of the k± jet-clustering algo- 
rithm for hadron colliders p0|-p2|. The main differences 
concern how particles are merged together and when the 
clustering stops. The different types of merging, or re- 
combination, schemes were investigated in Ref. pOf . D0 
chooses the scheme that corresponds to four-vector addi- 
tion of momenta, because pO[ : 

1. it is conceptually simple; 

2. it corresponds to the scheme used in the k±^ algo- 
rithm in e+e^ annihilations ||l^,|l^; 

3. it has no energy defect p2| , a measure of pertur- 
bative stability in the analysis of transverse energy 
density within jets; and 

4. it is better suited |3^] to the missing transverse en- 
ergy calculation in the jet-momentum calibration 
method used by D0. 

To stop clustering, D0 has adopted the proposal |2|] 
that halts clustering when all the jets are separated by 
AT?. > D. This rule is simple, and maintains a similarity 
with cone algorithms for hadronic collisions. The value 
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(a) 



•Beam 




2. If the minimuni of all possible da and dij is a dij, 
then replacing objects i and j by their merged object 



ij,p.y), where 



-Eij — Ei 




Pij = Pi + Pj- 

And if the minimum is a da, then removing object i from 
the list and defining it to be a jet. 

3. Repeating Steps 1 and 2 when there are any objects 
left in the list. 

The algorithm produces a list of jets, each separated 
by A7^ > D. Figure | illustrates how the fc_L algorithm 
successively merges the particles in a simplified diagram 
of a hadron collision. 



B. Preclustering 



FIG. 4. A simplified example of the final state of a col- 
lision between two hadrons. (a) The particles in the event 
(represented by arrows) comprise a list of objects, (b-f) Solid 
arrows represent the final jets reconstructed by the fcx algo- 
rithm, and open arrows represent objects not yet assigned to 
jets. The five diagrams show successive iterations of the al- 
gorithm. In each diagram, a jet is either defined (when it 
is well-separated from all other objects), or two objects are 
merged (when they have small relative fcx). The asterisk la- 
bels the relevant object(s) at each step. 



D = 1.0 treats initial-state radiation in the same way as 
final-state radiation |jll|,Q. 

The jet algorithm starts with a list of preclusters 
as defined in the next section. Initially, each preclus- 
ter is assigned a momentum four-vector {E, p) = 
^'preciuster(l,sin0cos(/!),sin0sin(^, COS0), written in terms 
of the precluster angles and The execution of the jet 
algorithm involves: 

1. Defining for each object i in the list: 



Px.i 



p 



and for each pair (i, j) of objects: 
A7^2. 

dij = mm [Pta,Pt.j\ 



r 2 2 

[PT,z,PT,j 



- Vj? + 



D2 



(3.2) 



where D is the stopping parameter of the jet algorithm. 
For D — l.Q and A7?.y <C 1, dij reduces to the square of 
the relative transverse momentum {k_\_) between objects. 



In the computer implementation of the jet algo- 
rithm, the processing time is proportional to , where 
A'^ is the number of particles (or energy signals) in the 
event |2^ . The zero-suppression circuit reduces the num- 
ber of calorimeter cells that have to be read out in each 
event. To reduce this further, we employ a preclustering 
algorithm. The procedure assigns calorimeter cells (or 
particles in a Monte Carlo event generator) to preclus- 
ters, suitable for input to the jet-clustering algorithm. In 
essence, calorimeter cells are collapsed into towers, and 
towers are merged if they are close together in (77, (j)) space 
or if they have small pr- Monte Carlo studies have shown 
that such preclustering reduces the impact of ambiguities 
due to calorimeter showering and finite segmentation, es- 
pecially on the reconstructed internal jet substructure. 
For example, when a single particle strikes the bound- 
ary between two calorimeter towers, it can produce two 
clusters of energy. Conversely, two collinear particles will 
often shower in a single calorimeter tower. In both cases, 
there is a potential discrepancy in the number of energy 
clusters found at the calorimeter level and the particle 
level. Preclustering at both the calorimeter and at the 
particle level within a radius larger than the calorimeter 
segmentation integrates over such discrepancies. 

The preclustering algorithm consists of the following 
six steps: 

1. Starting from a list of populated calorimeter cells 
in an event, remove any cells with Et < —0.5 GeV. Cells 
with such negative Et — rarely observed in minimum- 
biag|^ events (see Fig. |]) — are considered spurious. 



^The minimum- bias trigger requires a coincidence signal 
in the scintillating-tile hodoscopes located near the 

beampipe. 
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Energy (MeV) 



FIG. 5. Mean energies in calorimeter cells for a sample 
of minimum-bias events. The contribution from instrumen- 
tal effects is included, which occasionally leads to negative 
energy readings. For each cell, the energy distribution illus- 
trated in Fig. ti is fitted to a Gaussian. Before readout, the 
zero-suppression circuit in each cell's electronics sets to zero 
energy the channels in a symmetric window about the mean 
pedestal. These channels are not read out, causing the dip 
observed near zero. 



2. For each calorimeter cell centered at some {9, (f) rel- 
ative to the primary interaction vertex, define its pseu- 
dorapidity: 



77 



In tan ■ 



3. For each calorimeter tower sum the transverse 
energy of cells c in that tower: 

c et 

where Ec is the energy deposited in cell c. 

4. Starting at the extreme negative value of 77 and 
(/) = 0, combine any neighboring towers into preclus- 
ters such that no two preclusters are within A7?.p'° = 
y^Ari^ + A02 = 0.2. 
mass prescription p 



The combination follows the Snow- 



Ej" — Ef.i ~\' Ex.- 



V = 



Eta + Er.j 




0.15 0.20 

pre , 



E7(GeV) 



FIG. 6. The mean number of preclusters per event, as a 
function of the setting of minimum transverse energy required 
for preclusters [E'^'^). 

5. Because of pile-up in the calorimeter, precluster 
energies can fluctuate in both positive and negative di- 
rections. Preclusters that have negative transverse en- 
ergy Et = Et- < 0, are redistributed to k neighboring 
preclusters in the following way. Given a negative Et 
precluster with {Et- ,ri-,(j)-), we define a square S of 
size (77- ± 0.1) X {(f)- ± 0.1). When the following holds: 



J2etAv,<I>) > \Et- 



(3.3) 



k e S 



where only preclusters with positive Et that are located 
within the square S are included in the sum, then Et- 
is redistributed to the positive preclusters in the square, 
with each such precluster k absorbing a fraction 



E' 



T.k 



^Et,} 

k e S 



of the negative Et- If Eq. (3.3) is not satisfied, the 
"search square" is increased in steps of A77 — ±0.1 and 
Acjj = ±0.1, and another redistribution is attempted. 
If redistribution still fails for a square of {ri_ ±0.7) x 
{(f>- ± 0.7), the negative energy precluster is isolated in 
the calorimeter and ignored (by setting Et- = 0). 

6. Preclusters with < Et < E^'° = 0.2 GeV, are 
redistributed to neighboring preclusters, as specified in 
Step 5. To reduce the overall number of preclusters, we 
also require that the search square have at least three 
positive Et preclusters. The threshold EJ^'^^ was tuned 
to produce about 200 preclusters/event (see Fig. ^, in 
order to fit our constraints for processing time. 



Eta 



Et, 



3V] 



Eta + Et- 



C. Calibration of jet momentum 



The procedure evolves in the direction of increasing 0, 
and then increasing 77. 



A correct calibration of jet momentum reduces overall 
experimental uncertainties on jet production. The cali- 
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bration at D0 also accounts for the contribution of the 
underlying event (momentum transferred as a result of 
the soft interactions between the remnant partons of the 
proton and antiproton) . All such corrections enter in the 
relation between the momentum of a jet measured in the 
calorimeter p™'''*^ and the "true" jet momentum p''^"'' fs^j 



true 
Pjot 



JCt 



Poiv 



Jet 



(3.4) 



where po denotes an offset correction, i?jet is a correction 
for the response of the calorimeter to jets, and C is the 
instantaneous luminosity. A true jet is defined as being 
composed of only the final-state particle momenta from 
the hard parton-parton scatter (i .e., before interaction in 
the calorimeter). Although Eq. (3.4) is valid for any jet 
algorithm, po and the components of i?jet depend on the 
details of the jet algorithm. Our calibration procedure 
attempts to correct calorimeter-level jets (after interac- 
tions in the calorimeter) to their particle-level (before the 
individual particles interact in the calorimeter), using the 
described k± jet algorithm, with D — 1.0. The procedure 
follows closely that of calibration of the fixed-cone jet al- 
gorithm . The fixed-cone jet algorithm requires an 
additional scale factor in Eq. ( |3.4| ), but we find no need 
for that kind of calorimeter-showering correction in the 
k± jet momentum calibration psf . 

The offset po corresponds to the contribution to the 
momentum of a reconstructed jet that is not associated 
with the hard interaction. It contains two parts: 



PO = Ou 



O. 



where Oue is the offset due to the underlying event, 
and Ozb is an offset due to the overall detector envi- 
ronment. Ozb is attributed to any additional energy in 
the calorimeter cells of a jet from the combined effects of 
uranium noise, multiple interactions, and pile-up. The 
contributions of Oue and Ozb to k±^ jets are measured sep- 
arately, but using similar methods. The method overlays 
D0 data and Monte Carlo events, as described in what 
follows. 

The Monte Carlo events are generated by HERWIG 
(version 5.9) with 2^2 parton pr-thresholds of 
30, 50, 75, 100, and 150 GeV, and the underlying-event 
contribution switched off. The Monte Carlo events are 
propagated through a GEANT-based |3^ simulation of 
the D0 detector, which provides a cell-level simulation 
of the calorimeter response and resolution. These Monte 
Carlo events are then passed through the calorimeter- 
reconstruction and jet-finding packages, defining the ini- 
tial sample of jets. Detector simulation does not include 
the effects of uranium noise nor of the accelerator con- 
ditions causing multiple interactions and pile-up. The 
total contribution from these three effects is modeled us- 
ing zero-bias events, which correspond to observations 
at random pp bunch crossings. Zero-bias events were 
recorded by the D0 detector at different instantaneous 
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FIG. 7. The offset correction Ozb as a function 

of pseudorapidity of k±_ jet [D — 1.0). The offset 
Ozb accounts for the combined effects of pile-up, uranium 
noise, and muffiple interactions. The different sets of 
points are for events with different instantaneous luminos- 
ity C « 14, 10, 5, 3, 0.1 X 10^''cm"^s"\ The curves are fits to 
the points at different £, using the same functional form as 
employed for the cone algorithm in Ref. |3q]. 



luminosities in special data-taking runs without the zero- 
suppression discussed in Sec. The cell energies in zero- 
bias events are added cell-by-cell to the energies in sim- 
ulated Monte Carlo jet events. The summed cell ener- 
gies are then zero-suppressed offfine, using the pedestals 
appropriate to the zero-bias running conditions. Fi- 
nally, the summed cell energies are passed through the 
calorimeter-reconstruction and jet-finding packages, pro- 
ducing a second sample of jets. The two samples are 
compared on an event-by-event basis, associating the jets 
in events of the two samples that have their axes sepa- 
rated by ATZ < 0.5 jss). The difference in the measured 
Pt of the corresponding matched jets is Ozb, and shown 
in Fig. as a function of rj^^*", for different instantaneous 
luminosities. 

The event-overlay method was checked with the fixed- 
cone jet algorithm for TZ = 0.7. For jets with 30 GeV 
< Et < 50 GeV, this method gives only 14% (28%) 
smaller offsets [AOzb = 0.25 (0.39) GeV per jet], at 
£ « 5 (0.1) X lO^Ocm-^s-i relative to Ref. H. Indepen- 
dent of jet Et, the method used in Ref. measures the 
Et per unit A77 x A(j> in zero-bias events, and scales the 
value by the area of the jet cone. In the event-overlay 
method, Ozb decreases by as much as 40% when the 
cone-jet transverse energy increases to 125 GeV < Et < 
170 GeV. Approximately 30% of this decrease can be ex- 
plained by the i?;^*-dependence of the occupancy of cells 
within cone jets (the fraction of cells with significant en- 
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ergy deposition inside the cone). The remaining 70% of 
the Ozh dependence on jet Et is assigned as a systematic 
uncertainty on our method. Since the observed depen- 
dence is less pronounced in the k± jet algorithm, this 
error amounts at most to 15% in the highest jet pt bin. 
In addition, wc include a systematic uncertainty of 0.2 
GeV arising from the fits in Fig. |^. Using our overlay 
method for both algorithms, the offsets Ozh in the k± 
jet algorithm (with D = 1.0) are generally 50 — 75% (or 
about 1 GeV per jet) larger than in the fixed-cone jet 
algorithm (with 7^ = 0.7) (3|]. 

The offset due to the underlying event Ouc is modeled 
with minimum-bias events. A minimum-bias event is a 
zero-bias event with the additional requirement of a co- 
incidence signal in the scintillating-tile hodoscopes ||l9| ] 
near the beampipe. The additional requirement means 
there was an inelastic pp collision during the bunch cross- 
ing. In addition to Ouc, a minimum-bias event in the 
D0 calorimeter includes energy from uranium noise, mul- 
tiple interactions, and pile-up. However, the luminos- 
ity dependence of multiple interactions and pile-up in 
minimum-bias events is different than in zero-bias events. 
In the limit of very small luminosity, these contributions 
are negligible, and a minimum-bias event at low lumi- 
nosity therefore contains the offset due to the underly- 
ing event and uranium noise, while a zero-bias event at 
low luminosity has only the offset from uranium noise. 
To measure Ouc, we again compare two samples ofjets. 
Minimum-bias events as measured by the D0 calorimeter 
at low luminosity are added to Monte Carlo jet events, 
where the resulting jets define the first sample of jets 
in the determination of Oue- The second sample ofjets 
is reconstructed from zero-bias events at low luminosity 
and also added to Monte Carlo jet events. On an event- 
by-event basis, Ouc is calculated by subtracting the mo- 
mentum of jets in the second sample from the momen- 
tum of matching jets in the first sample. The underlying 
event offset Ouc for k± jets is shown in Fig. ||. Using this 
method for both algorithms, the offset Ouc for k± jets 
(with D = 1.0) is found to be approximately 30% larger 
than for the fixed-cone jet algorithm (with TZ = 0.7). 

D0 measures the jet momentum response based on 
conservation of pt in photon-jet (7-jet) events [Q. The 
electromagnetic energy /momentum scale is determined 
from the Z,J/ip — > e~''e~, and ir'^ — > 77 — s- e+e^e^e^ 
data samples, using the known masses of these particles. 
For the case of a 7-jet two-body process, the jet momen- 
tum response can be characterized as: 



R 



'jet 



(3.5) 



where pxj and n are the transverse momentum and di- 
rection of the photon, and ]^rp is the missing transverse 
energy, defined as the negative of the vector sum of the 
transverse energies of the cells in the calorimeter. To 
avoid resolution and trigger biases, i?jet is binned in terms 
of E' = p^;^*^'*-cosh(r;jet) and then mapped onto the offset- 
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FIG. 8. The correction for underlying event Oue as a 
function of |j7| for k_i_ jets {D — 1.0). The solid curve is the 
fit of the results for the cone jet algorithm in Ref. |^ scaled 
to the results for the k± jet algorithm. The dashed curves 
denote the one standard deviation (s.d.) systematic error. 



corrected jet momentum. Thus, E' depends only on pho- 
ton variables and jet pseudorapidity, which are quantities 
that are measured with very good resolution. i?jct and 
E' depend only on the jet position, which has little de- 
pendence on the type of jet algorithm employed. 

i?jot as a function of p'^°'^^ — Po for kj_ jets is shown 
in Fig. 0. The data points are fitted with the functional 



form i?jct(p) 



b\n{p) 



nn(p))^. The response i?jct 
for cone jets (with TZ — 0.7) ||3q] and for k± jets {D — 1.0) 
is different by about 0.05. This difference does not have 
any physical meaning; it corresponds to different voltage- 
to-energy conversion factors at the cell level used in the 
reconstruction of jets. 



D. Comparison of the k± jet algorithm to the cone 
jet algorithm 

It is of interest to compare the momenta of k± jets 
to those ofjets reconstructed with the D0 fixed-cone al- 
gorithm |30|. These results refer to the k± jet algorithm 
described above with D = 1.0 and co rrected according to 
the prescriptions given in Sec. Ill C . The cone jets were 
reconstructed with TZ = 0.7 and corrected according 
to Ref. 1^^. This comparison involves about 75% of the 
events in the 1994-1996 data that were used for the anal- 
ysis of the inclusive cone-jet cross section at y/s — 1800 
GeV [P3 ■ The two algorithms are similar by design [E2| , 
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FIG. 9. The response correction for k± jets with D — 1.0, 
as a function of offset-corrected jet momentum. The Monte 
Carlo point (*) is used to constrain the fit (solid) at high pj^t'"'- 
The dashed curves denote the ±1 s.d. systematic error. 

defining similar jet directions and momenta, at least for 
the two leading (highest pt) jets in the event. The re- 
maining jets in the event usually have much smaller px, 
making them more difficult to measure, and so we do 
not consider them here. The jets reconstructed by each 
algorithm are compared on an event-by-event basis, as- 
sociating a cone jet with a k± jet if they are separated 
by An < 0.5. 

To obtain a sample of events with only good hadronic 
jets, the following requirements were placed on the events 
and on the leading two reconstructed k± jets. These 
criteria are based on standard jet quality requirements 
(to remove spurious clusters) in use at D0 for the fixed- 
cone jet algorithm [p7| : 

• Measured event vertex was required to be within 
50 cm of the center of the detector. 

• l^rpl was required to be less than 70% of the pr of 
the leading jet. 

• Fraction of jet px measured in the coarse hadronic 
calorimetry was required to be less than 40% of the 
total jet Pt- 

• Fraction of jet px measured in the electromagnetic 
calorimetry was required to be between 5% and 
95% of the total jet px- 

• Jets were required to have \ri\ < 0.5. 

These requirements yield a sample of 68946 k± jets. 
The axes of 99.94% of these jets are reconstructed within 



FIG. 10. The distance A7^ = ^/Arf+Atj^ between a 
fcx-jet axis and its matching cone-jet axis. The k± jets were 
reconstructed with D = 1.0, and the cone jets were recon- 
structed with TZ = 0.7. Only the two leading jets from each 
algorithm were considered. The k± jets were selected with 
\t]\ < 0.5. 



ATI < 0.5 of a cone-jet axis, when the matching jet 
is one of the two leading cone jets in the event. For 
such pairs of jets, the distance between a fcjjet axis 
and matching cone-jet axis is shown in Fig. |l^. The 
fixed-cone algorithm finds a jet within ATZ < 0.1 of a 
k± jet 91% of the time. Figure ^ shows the differ- 
ence pt(/c_l jet) — i?T (cone jet) as a function of pT(fc_L jet). 
Generally, the px of k± jets (D = 1.0) is higher than the 
Ex of associated cone jets {TZ = 0.7). The difference in- 
creases approximately linearly with jet px, from about 5 
GeV (or 6%) at px « 90 GeV to about 8 GeV (or 3%) at 
Px ~ 240 GeV. This may be explained by how the two 
algorithms deal with hadronization effects p^. 



E. Subjets 

The subjet multiplicity is a natural observable for char- 
acterizing a k± jet 21|. Subjets are de fined by reap- 
plying the k± algorithm, as in Sec. Ill A, starting with 
a list of preclusters assigned to a particular jet. Pairs 
of objects with the smallest dij are merged successively 
until all remaining pairs of objects have 



2 2 ^^?7 2 

dij = mm{px^i,Pxj)^j2~ > VcutPxH'^^): 



(3.6) 
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FIG. 11. The difference pT{k± jet) — Ericone jet) as a 
function of the k± jet pr- A cone jet is associated with a k± 
jet if their axes are separated by ATZ < 0.5. 



where PtO^I) is the px of the entire jet in the k± algo- 
rithm described above, and < ?/cut < 1 is a dimension- 
less parameter. Objects satisfying Eq. (3.6) are called 



subjets, and the number of subjets is the subjet multi- 
plicity M of a jet. For ycut = 1, the entire jet con- 
sists of a single subjet (M = 1). As j/cut decreases, the 
subjet multiplicity increases, until every precluster be- 
comes resolved as a separate subjet in the limit ?/cut — > 0. 
Two subjets in a jet can be resolved when they are not 
collinear (i.e., well-separated in rj x (j) space), or if they 
are both hard (i.e., carry a significant fraction of the jet 
Pt)- 

We now turn to the theoretical treatment of sub- 
jet multiplicity. Perturbative and resummed calcula- 



tions [ |liyi7| and Monte Carlo estimates (see Sec. IVD) 
predict that gluon jets have a higher mean subjet multi- 
plicity than quark jets. To understand the origin of this 
prediction, we consider first how a jet can contain multi- 
ple subjets. Clearly, at leading-order, 2 — 2 subprocesses 
yield M — 1. However, higher-order QCD radiation can 
increase the average value of M. At next-to-leading or- 
der, there can be three partons in the final state of a pp 
collision. If two partons are clustered together into a jet, 
they can be resolved as distinct subjets (M — 2) for a 
sufficiently small choice of j/cut- For larger ?/cut, the value 
of M depends on the magnitude and direction of the ra- 
diated third parton. In QCD, the radiation of a parton 
is governed by the DGLAP splitting functions [^|. The 
radiated third parton is usually soft and/or collinear with 
one of the other two partons, leading to jets with Af = 1. 



However, hard or large-angle radiation, although rare, 
causes some jets to have M = 2. Consequently, when 
many jets are analyzed using some high j/cut, the two- 
subjet rate will yield (M) > 1. 

In the framework of parton showers, repeated appli- 
cation of DGLAP splitting provides jets with M > 2. 
Monte Carlo event generators incorporate parton show- 
ers into the initial and final states of a 2 ^ 2 hard scatter. 
Because of its larger color factor, a parton shower initi- 
ated by a gluon in the final state will tend to produce a 
jet with more subjets than one initiated by a quark. Sim- 
ilarly, a soft parton radiated in the initial state will tend 
to cluster with a hard final-state parton when ATZ < D. 
For the case of initial-state radiation, the subjet multi- 
plicity depends weakly on whether the final-state par- 
tons in the 2 — > 2 hard scatter are quarks or gluons. 
The contribution of initial-state radiation to the subjet 
multiplicity does, however, depend on y/s. Initial-state 
radiation is treated on an equal footing as final-state ra- 
diation in the k± algorithm with D — 1.0 |l^,Q, and 
diminishes in importance as D decreases. In general, 
subsequent emissions in parton showers have less energy 
and momentum, and this structure is revealed at smaller 
2/cut values through an increase in the subjet multiplicity: 
(A^(2/cut)) > (Af(2/cut)), where y'^^, < T/cut. 

Experimentally, the growth of M at very small j/cut is 
reduced by the granularity of the detector and by the 
preclustering algorithm. Theoretical predictions for M 
are therefore treated in the same way as the experimen- 
tal measurements, i.e., by preclustering (as in Sec. HI E| ). 
Requiring preclusters to be separated by AT^p'''', means 
that the subjets nearest in (77, (j)) space begin to be re- 
solved for 



ycut < 



/A7^P'■ 



V 2D 



(3.7) 



based solely on the fraction of pr carried by the sub- 
jet in the jet. The factor 1/2 corresponds to the maxi- 
mum fraction of jet Pt carried by the softest subjet [see 
Eq. (3.6)]. The preclustering stage provides a comparison 
of the measurement of M with prediction in the interest- 
ing region of small ?/cut , without an explicit correction for 
detector granularity. 

The subjet analysis in this paper uses a single resolu- 
tion parameter j/cut = 10""^. For this ?/cut, the minimum 
subjet Pt is approximately 3% of the total jet pT, in- 
dependent of the choice of the D par ameter. Because 
?/cutj as defined by Eqs. (3.2) and (3.6), involves a ratio 
of subjet pt to jet px, the subjet multiplicity is there- 
fore not significantly sensitive to multiplicative changes 
in the overall pt scale. Consequently, given the fact that 
subjets are specified during jet reconstruction, and the jet 
momentum calibration is derived after reconstruction, we 
do not attempt to correct the momenta of individual sub- 
jets. However, the subjet multiplicity is corrected for the 
experimental effects that cause an offset in jet px- In gen- 
eral, the presence of uranium noise, multiple interactions, 
and pile-up, tends to increases the subjet multiplicity. 
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IV. DATA SAMPLES 

In leading-order QCD, the fraction of final-state jets 
originating from gluons decreases with increasing x oc 
Pt/V^, the momentum fraction carried by the initial- 
state partons. This is due primarily to the x-dependence 
of the parton distributions. Because, for fixed px, the 
gluon fraction decreases when ^/s is decreased from 1800 
GeV to 630 GeV, this suggests an experimental way to 
define jet samples with different mixtures of quarks and 
gluons. A single set of criteria can be used to select jets at 
the two beam energies, without changing any of the de- 
tector elements. We use this principle to analyze an event 
sample recorded at the end of 1995 by the D0 detector at 
= 630 GeV, and compare it with the larger 1994-1995 
event sample collected at i/i = 1800 GeV. The lower 
range of jet px populated by the smaller event sample 
at y/s = 630 GeV dictat ed th e ultimate criteria used in 
the comparison. In Sec. IVA, we first describe a simple 



test of a set of criteria used to sele ct qua rk-enriched and 
gluon-enriched jet samples. In Sec. [V B|, w e specify each 
criterion used in the analysis. In Sec. |[V C , we provide a 
Monte Carlo estimate of the quark/gl uon yi eld based on 
the full set of criteria. Finally, in Sec. [VP , we describe 
how to estimate the subjet content of gluon and quark 
jets. 



A. Gluon and quark samples at leading-order in 
QCD 

For a given set of parton distribution functions (PDFs), 
the relative admixture of gluon and quark jets passing a 
set of kinematic criteria can be estimated using a leading- 
order QCD event generator. At this order, there is no 
dependence on jet algorithm, because each of the two 
final-state partons defines a jet. We use the HERWIG 
v5.9 Monte Carlo with the CTEQ4M || PDFs to gen- 
erate leading-order QCD 2^2 events, and keep track of 
the identity (gluon or quark) of the partons. At leading 
order, the gluon-jet fraction / corresponds to the num- 
ber of final-state gluons that pass the selections divided 
by the total number of final-state partons that pass the 
selections. For example, the jet sample selected from 
only gg ^ gg or qq — > gg events will have a gluon-jet 
fraction of unity. Figure O shows that for the full en- 
semble of Monte Carlo events, the gluon-jet fraction at 
y/s = 630 GeV is about 30% smaller than at ^ = 1800 
GeV, where we have selected central {\ri\ < 0.5) jets with 
minimum parton px ~ 55 GeV and maximum parton 
Pt = 100 GeV. This difference is due primarily to the 
relative abundance of initial-state gluons at these x val- 
ues for y/s — 1800 GeV compared to y/s = 630 GeV. 



H- 0.8 



o 

CC 0.7 



CD 



O 

O 



0.6 



0.4 



0.3 



• 1800 GeV 
o 630 GeV 




52 54 56 58 

Minimum parton p-p (GeV) 



FIG. 12. The Monte Carlo gluon-jet fraction / at lead- 
ing-order, for final-state partons with maximum parton 
Pt = 100 GeV, and minimum parton pr ~ 55 GeV, as a func- 
tion of the minimum parton pr, using the CTEQ4M PDF. 
Both partons are required to be central {\ri\ < 0.5). The solid 
symbols show the prediction for ^/s = 1800 GeV, and the 
open symbols show the prediction for -^s — 630 CeV. 



B. Jet data samples 

We define gluon-enriched and quark-enriched central 
(|?7| < 0.5) jet samples using identical criteria at — 
1800 GeV and 630 GeV, thereby reducing any experi- 
mental biases and systematic effects. We select events 
that pass a trigger requiring the scalar sum of above 
30 GeV within a cone of si ze TZ = 0.7 and apply 

the selections listed in Sec. HID, but only for jets with 



measured pt between 55 and 100 GeV. These cuts yield 
samples of 11,007 jets at ^/i = 1800 GeV, and 1194 jets 
at Vs = 630 GeV. 

An important point is that these jets were recon- 
structed with the k± algorithm for D — 0.5. This choice 
tends to select events with fewer subjets from initial-state 
radiation, which can vary with ^/s (see Sec. HIE). Fig- 
ure |l^ shows that the pt distribution of the selected jets 
at = 1800 GeV is harder than at = 630 GeV. 
The mean jet pr at y/^ = 1800 GeV is 66.3 ±0.1 GeV, 
which is 2.3 GeV higher than at y/s = 630 GeV. This 
cannot be caused by any differences in the contribution 
to the offset in the jet pr- In fact, the entire offset is 
po « 3 - 4 Ge V per jet at = 1800 GeV for D = 1.0 
(see Sec. Ill C ), and is therefore an expected factor 4 
smaller for D = 0.5. Moreover, only a small fraction of 
the jet offset can be attributed to the difference in y/s. 
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FIG. 13. ThepT distribution of selected central {\ri\ < 0.5) 
jets in D0 data, before applying a cutoff on jet pr- The data 
at = 630 GeV are normalized to the data at = 1800 
GeV in the bin 54 < pt < 60 GeV. The turnover at lower 
jet pt is due to inefficiencies in the trigger. For the following 
analysis, we use jets with 55 < pt < 100 GeV. 



Even so, offset differences can only change the sub jet 
multiplicity by shifting the relative jet px- Rather than 
attempting to measure and account for such small effects 
in the jet pt distributions, we simply use identical jet 
criteria at the two beam energies, and estimate the un- 
certainty on AI by varying the jet selection cutoffs (see 
Sec. |VC|). 



C. Jet samples in Monte Carlo events 

To estimate the number of gluon jets in the y/s — 1800 
GeV and 630 GeV jet samples, we generated approxi- 
mately 10,000 HERWIG events at each y/s, with par- 
ton Pt > 50 GeV, and requiring at least one of the 
two leading-order partons to be central {\ri\ < 0.9). The 
events were passed through a full simulation of the D0 
detector. To simulate the effects of uranium noise, pile- 
up from previous bunch crossings, and multiple pp in- 
teractions in the same bunch crossing, we overlaid D0 
random-crossing events onto our Monte Garlo sample, 
on a cell-by-cell basis in the calorimeter. (A sample with 
instantaneous luminosity of £ « 5 x lO'^^cm^^s^^ was 
used at y/i = 1800 GeV, and £ « 0.1 x lO^^cm-^s"! was 
used at -^/s — 630 GeV.) These pseudo events were then 
passed through the normal offline-reconstruction and jet- 
finding packages. Jets were then selected using the same 
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FIG. 14. The normalized pt distribution of central 
{\Tj\ < 0.5) jets selected in Monte Carlo events at -^s = 1800 
GeV and 630 GeV. Each distribution has been normalized to 
unit area. 



criteria as used for D0 data, and their px distribution is 
shown in Fig. |lj. 

We tag each such selected Monte Carlo jet as either 
quark or gluon based on the identity of the nearer (in 
1] X (p space) final-state parton in the QCD 2^2 hard 
scatter. The distance between one of the partons and 
the closest calorimeter jet is shown in Fig. |l^. There 
is clear correlation between jets in the calorimeter and 
partons from the hard scatter. The fraction of gluon 
jets is shown in Fig. |l^ as a function of the minimum 
Pt used to select the jets. There is good agreement for 
the gluon-jet fraction obtained using jets reconstructed 
at the calorimeter and at the particle levels (A/ < 0.03). 
The smaller gluon-jet fractions relative to leading-order 
(Fig. |lj) are due mainly to the presence of higher-order 
radiation in the QGD Monte Carlo. When pT cutoffs 
are applied to particle-level jets, the associated leading- 
order partons shift to significantly higher pT- Since the 
gluon-jet fraction decreases with increasing parton pT, f 
is smaller when events are selected according to particle- 
level jet Pt rather than when they are selected according 
to partonic pt- The same is true for cutoffs applied to 
the calorimeter-level jets compared to the particle-level 
jets, although here the A/ discrepancy is much smaller. 
In what follows, we shall use nominal gluon-jet fractions 
/i800 — 0.59 and /gso = 0.33, obtained from Monte Carlo 
at the calorimeter level for 55 < pr < 100 GeV. 
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FIG. 15. The distance of the closest calorimeter-level 
Monte Carlo jet to one of the leading final-state partons. 
The solid (open) points show the Monte Carlo sample at 
= 1800 (630) GeV. Each distribution has been normal- 
ized to unit area. 
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D. Sub jets in gluon and quark jets 

Using the previously described jet samples, there is a 
simple way to distinguish between gluon and quark jets 
on a statistical basis |2^. The subjet multiplicity in a 
mixed sample of gluon and quark jets can be written as 
a linear combination of subjet multiplicity in gluon Mg 
and quark jets Mq-. 



M = fMg + (1 - f)Mq 



(4.1) 



The coefficients are the fractions of gluon and quark jets 
in the mixed sample, / and (1 — /), respectively. Con- 
sidering Eq. ([4.1|) for two samples of jets at ^/s = 1800 



GeV and 630 GeV, and assuming that Mg and Mq are 
independent of (we address this assumption later), 
we can write: 



(1 — /eao) Afisoo 



/l80o) -^'^630 



/: 



1800 



7630 



Mq = 



/l800-^^630 



/630M: 



1800 



/: 



1800 



630 



(4.2) 



(4.3) 



where Afiggo and Mq^q are the measured multiplicities in 
the mixed-jet samples at y/s — 1800 GeV and 630 GeV, 
and /i8oo and /gao are the gluon-jet fractions in the two 
samples. The extraction of Alg and Mq requires prior 
kno wledge of the two gluon-jet fractions, as described in 



Sec. [VC, Since the gluon-jet fractions depend on jet 
Pt and 77, Eqs. ( |4.2D and ( |4.3| ) hold only within restricted 
regions of pha se sp ace, i .e., over small ranges of jet pt and 
77. Equations (4^) and (O) can, of course, be generalized 



to any observable associated with a jet. 

We use our Monte Carlo samples to check Eqs. (|4.2| ) 
and (^^) for k± jets reconstructed using the full-detector 
simulation with D = 0.5. Such a consistency test does 
not depend on the details of the subjet multiplicity distri- 
butions {Mq, Mg, MisQo, Mq3q). The extracted distribu- 
tions in Mg and Mq are shown in Fig. |l^. As expected, 
Monte Carlo gluon jets have more subjets, on average, 
than Monte Carlo quark jets: (Mg) > (Mq). This is also 
found for jets reconstructed at the particle level, and the 
differences between gluon and quark jets do not appear 
to be affected by the detector. Also, the subjet multi- 
plicity distributions for tagged jets are similar at the two 
center-of-mass en ergi es, ver ifyin g the assumptions used 
in deriving Eqs. (^) and (|4.3| ). Finally, the extracted 
Alq and A'lg distributions agree very well with the tagged 
distributions. This demonstrates self-consistency of the 
extraction using Eqs. (|4.2D and (|4.3|). 



FIG. 16. The gluon-jet fraction of selected jets with max- 
imum Pt ~ 100 GeV and minimum pr between 52 and 58 
GeV, as a function of minimum jet pt, for y's = 1800 GeV 
and 630 GeV, using the CTEQ4M PDF. The jets have been 
tagged through the identity of the nearer leading-order fi- 
nal-state parton. 



V. SUBJET MULTIPLICITIES 

A. Uncorrected subjet multiplicity 

Figure ^ shows the distributions of subjet multiplicity 
for the D0 data samples described in Sec. This is the 
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FIG. 17. Uncorrected subjet multiplicity in fully-simulated 
Monte Carlo (a) gluon and (b) quark jets. The number of 
jets A'jots(Af) in each bin of subjet multiplicity on the ver- 
tical axis is normalized to the total number of jets in each 
sample A'j';,"^ = X^a/ -^j<=ts(-^)- The measured distributions 
(solid) are extracted from the mixed Monte Carlo jet sam- 
ples at = 1800 GeV and 630 GeV. The tagged distribu- 
tions (open) are for -^s = 1800 GeV (triangles) and 630 GeV 
(squares). 



first measurement of its kind at a hadron collider. The 
average number of subjets in jets at -y/s = 1800 GeV 
is (Afisoo) = 2.74 ± 0.01, where the error is statistical. 
This is higher than the value of (Meao) = 2.54 ± 0.03 at 
\fs = 630 GeV. The observed shift is consistent with the 
prediction that there are more gluon jets in the sample 
at y/s = 1800 GeV than in the sample at ^/s = 630 
GeV, and that gluons radiate more subjets than quarks 
do. The fact that the spectrum is harder at yfs = 
1800 GeV than at y/s — 630 GeV cannot be the cause 
of this effect because the subjet multiplicity decreases 
with increasing jet pt- Figure ^ shows the rather mild 
dependence of the average subjet multiplicity on jet -pT- 

Subjets were defined through the product of their frac- 
tional jet pt a nd their separation in (r;, 0) space [see 
Eqs. (Ij) and (jsj)]. As shown in Figs. |2g and |l[ the 
shapes of the subjet spectra of the selected jets are 
similar at the two beam energies. The distributions sug- 
gest that jets are composed of a hard component and a 
soft component. The peak at about 55 GeV and fall-off 
at higher is due to single-subjet jets and the jet 
selections (55 < pr < 100 GeV). The threshold at subjet 
Pt ~ 1-75 GeV is set by the value ycut = 10^'^ and the 
minimum jet pt in the sample. 

While the Miggo and Mgao inclusive measurements 



FIG. 18. Uncorrected subjet multiplicity in jets from D0 
data at ^ = 1800 GeV and 630 GeV. 




80 90 100 

Jet p-p (GeV) 



FIG. 19. Uncorrected mean subjet multiplicity versus jet 
Pt in D0 data at ^/s — 1800 GeV. Note the suppressed zero 
on the vertical axis. 



at y/s — 1800 GeV and ^/s = 630 GeV are interest- 
ing in themselves, they can be interpreted in ter ms o f 
their glu on and quark content. According to Eqs. (4.2) 
and (43) the distributions in Fig. |l8| and their gluon- 
jet fractions at the two beam energies can yield the un- 
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FIG. 20. The uncorrected px distribution of subjets in 
data for jets with 55 < pr < 100 GeV and I?;! < 0.5. All 
selections have been applied, and each distribution has been 
normalized to unit area. 



FIG. 21. Same as in Fig. but with the low pt region 
expanded. The increase at low pr is observed for all j/cut, but 
the specific cutoff at pT(subjet) ~ 1.75 GeV is determined by 
our chosen value of j/cut = 10"'^. 



corrected subjet multiplicity distributions in gluon and 
quark jets. The extracted measurements of Mg and Mq 
are shown in Fig. ^ for the nominal values /isoo = 0.59 
and /eso = 0.33. As in the Monte Carlo simulation, the 
D0 data clearly indicate the presence of more subjets in 
gluon jets than in quark jets. Such distributions can be 
used directly (without correcting the subjet multiplici- 
ties) to discriminate between gluon and quark jets. The 
results depend only on Monte Carlo estimates of gluon- 
jet fractions at the two values of ^/s, and not on any 
detailed simulation of jet structure. 

The sensitivity of Mg and Mg to the assumed values of 
/i8oo and /gso was checked by investigating how the sig- 
nal (i.e., the difference between Mg and Mg) depended 
on this choice. It was found that when the gluon-jet 
fractions are either both increased or both decreased, 
the signal remains relatively unchanged. However, when 
the gluon-jet fractions are changed in opposite directions, 
this produces the largest change in the difference between 
gluon and quark jets. The result of using /isoo = 0.61 
and /eao = 0.30, instead of their nominal values, is shown 
in the extracted distributions of Fig. The Mg and Mq 
distributions of Fig. ^ are qualitatively similar to those 
of Fig. and the large difference between gluon and 
quark jets is still apparent. 

The subjet multiplicity distributions can be character- 
ized by their means (M), and by (M) — 1, which cor- 
respond to the average number of subjet emissions in a 
gluon or quark jet. For the nominal uncorrected D0 data 
shown in Fig. |2|, (M™'=^") = 3.05 ± 0.06 and (Af^"^''^") = 



2.28 ± 0.08. The analogous values for the Monte Carlo 
events (see Fig. ^ are (M™'=^") = 3.01 ± 0.09 and 
^^,jmeas^ =2.28±0.08. Becausc the quotcd Statistical Un- 
certainty on (M™°'''')is correlated with that on {M'^°'^), 
we define a ratio p3| , p^ of emissions in gluon jets to quark 
jets: 



{Mg 



{Mq 



(5.1) 



A value of r = 1 would mean that the substructure of 
gluon jets does not differ from that of quark jets. The ra- 
tio has a value of r = 1.61±0.15 for the uncorrected data 
of Fig. 11^, and r = 1.58 ± 0.16 for the analogous Monte 
Carlo events of Fig. |l^, where both uncertainties are sta- 
tistical. Using different values for gluon-jet fraction at 
the two values of ^/s (as in Fig. ^3|), yields the range of r 
values given in Table |. As expected, the observed ratio 
is smallest when the fraction of gluon jets increases at 
Vs = 1800 GeV and decreases at = 630 GeV. The 
two values of / are the only assumptions from Monte 
Carlo, and correspond to the largest source of system- 
atic uncertainty on r (described more fully in Sec. VC). 
In all cases, we find that r is significantly greater than 
unity, meaning that gluon jets and quark jets differ in 
their substructure. 
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FIG. 22. Uncorrected subjet multiplicity in gluon and 
quark jets, extracted from D0 data at = 1800 GeV and 
630 GeV, using nominal gluon-jet fractions /isoo = 0.59 and 
/630 = 0.33. 



/l800 


/630 
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0.59 


0.33 


3.05 ± 0.06 


2.28 ± 0.08 


1.61 ± 0.15 


0.61 


0.30 


2.99 ± 0.05 


2.34 ± 0.07 


1.49 ± 0.11 


0.61 


0.36 


3.05 ± 0.06 


2.24 ± 0.09 


1.65 ± 0.16 


0.57 


0.30 


3.06 ± 0.06 


2.31 ± 0.07 


1.57 ± 0.14 


0.57 


0.36 


3.15 ± 0.08 


2.19 ± 0.10 


1.81 ± 0.22 



TABLE I. The uncorrected subjet multiplicity in gluon and 
quark jets, and their ratio, extracted from D0 data, assuming 
different values of gluon-jet fractions at the two center-of-mass 
energies, based, in part, on Figs. |l^ and ^ 



B. Corrected subjet multiplicity 

As was stat ed ab ove, the experimental conditions de- 
scribed in Sec. p^II C smear the measurement of the subjet 
multiplicity. Although r expresses differences between 
gluon and quark jets as a ratio of mean subjet multi- 
phcities, the extracted Mg and Mq distributions need 



separate corrections for the various detector-dependent 
effects that can affect the value of r. The corrections are 
derived using Monte Carlo events, which are in agreement 
with the uncorrected D0 data, as shown in Figs. ^ and 
E3. The decomposition of the Monte Carlo events into 



Mg and Mq components was discussed in Sec. IV D. The 



distributions shown in Fig. ^ represent the uncorrected 
results for Monte Carlo events that we use to derive the 
unsmearing corrections. 

The corrected distributions of Alg and Mq are defined 
in Monte Carlo jets at the particle level (i.e., before devel- 



FIG. 23. Uncorrected subjet multiplicity in gluon and 
quark jets, extracted from D0 data at -^s — 1800 GeV 
and 630 GeV, using gluon-jet fractions /isoo = 0.61 and 
/630 = 0.30. 
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FIG. 24. Uncorrected subjet multiplicity in jets from D0 
and fully-simulated Monte Carlo events at ^/s = 1800 GeV. 



opment in the calorimeter). All selected calorimeter-level 
jets are matched (within ATI < 0.5) to jets reconstructed 
at the particle level. The matching procedure implicitly 
accounts for any mismeasurement of jet pt because there 
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FIG. 25. Uncorrected subjet multiplicity in jets from D0 
and from fully-simulated Monte Carlo events at ^/s — 630 
GeV. 



is no pt requirement in the matching. The preclustering 
and clustering algorithms applied at the particle level are 
identical to those applied at the detector level. We tag 
simulated detector jets as either gluons or quarks, and 
correlate the subjet multiplicity in particle jets (M*''"°) 
with that of detector partners {M'^°^^). These correla- 
tions are shown in Fig. at ^/s = 1800 GeV, and define 
the correction applied to the subjet multiplicity. Similar 
results are available at y/s = 630 GeV (not shown). 

The correction retrieves Af*'"° from M"^'^'^, in bins of 
^^meas general, the distributions of M*™<= and M*™'= 
in Fig. are shifted to lower values relative to M^'^'^^ 
and M™°'*'^. The shift in M is due mainly to the ef- 
fects of showering in the calorimeter, rather than from 
the combined effects of multiple interactions, pile-up, and 
uranium noise, which are reduced by using D = 0.5. For- 
tunately, shower development is independent of beam en- 
ergy, and the other contributions differ only slightly (see 
Sec. |VC| ). 

Shower development in the calorimeter tends to add 
subjets to a jet because any single particle can deposit 
energy in several towers of the calorimeter. Signals in 
many towers generate a large number of preclusters, and 
in turn, a large number of subjets. However, the oppo- 
site can also occur. For example, when two subjets at 
the particle level (each composed of one or two hadrons) 
deposit energy in a region of the calorimeter between 
them, such energy can "bridge" distinct subjets at the 
particle level into a single subjet at the calorimeter level. 
This bridging effect is more pronounced in jets that al- 



FIG. 26. The subjet multiplicity at particle-level (M*''"") 
versus the subjet multiplicity at calorimeter-level (M™°''°) (in- 
cludes effects of luminosity), at ^/s = 1800 GeV, for (a) gluon 
and (b) quark jets. 



ready have a large M*™". For this reason, the effects of 
multiple interactions, pile-up, and uranium noise tend to 
reduce the correction to M^'^^^ . 

To check that the correction defined by the correlations 
in Fig. is valid, it was applied to the uncorrected Mg 
and Mq Monte Carlo distributions in Fig. |l^. The re- 
sulting corrected distributions for Mg and Mq are given 
in Fig. and Fig. 28, respectively. The correction re- 



duces the average subjet multiplicity in the Monte Carlo 
to (M^™'=) = 2.19 ± 0.04 and (M^™<=) = 1.66 ± 0.04 and 
the corrected ratio is r = 1.82 ± 0.16. Any remaining 
small differences between the extracted and the tagged 
jy^truc distributions in Fig. ^ and Fig. |2^ are attributable 
to the differences between the extracted and the tagged 
^jmeas ^ ^ ^^ggO GcV) of Fig. |l|. Thcsc differences 
are smaller for the corrected distributions A/''"'', than 
for the uncorrected distributions. 

Figure |2^ shows the corrected subjet multiplicities for 
gluon and quark jets. The rate for M = \ quark jets 
has almost doubled, while the rate for M — 3 quark 
jets has fallen by a factor of ~ 2, relative to the un- 
corrected result. A similar effect is observed for gluon 
jets. From Fig. |2^, we obtain the corrected mean val- 
ues in the D0 data to be (Af*''™) = 2.21 ± 0.03 and 
(Aft'uc) = 1.69±0.04, which gives r = 1.75±0.15, in good 
agreement with the prediction from HERWIG. The un- 
smearing therefore widens the difference between gluon 
and quark jets. 

We choose not to correct M for any impact of the 
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FIG. 27. The subjet multiplicity in Monte Carlo gluon 
jets. The extracted distribution has been unsmeared. The 
tagged distribution was obtained directly from particle-level 
gluon jets at ^ = 1800 GeV. 



preclustering algorithm on subjet multiplicity. Instead, 
the preclustering algorithm can be applied easily to the 
particle-level events in Monte Carlo, and these are there- 
fore treated in the same way as the D0 data. For com- 
pleteness, we note that r can decrease by as much as 0.2 
at the particle level, when preclustering is turned off. 



C. Additional corrections and systematic 
uncertainties 

The dominant systematic uncertainty on the subjet 
multiplicity arises from the uncertainty on the gluon-jet 
fractions. In fixed-order perturbative QCD, the jet cross 
section at any given px is a more-steeply-falling function 
of pT at Vs = 630 GeV than at = 1800 GcV p?)- 
Consequently, applying identical cutoffs biases the {pr) 
of jets at ^/s = 1800 GeV upwards relative to ^/s ~ 630 
GeV. Monte Carlo studies indicate this bias is approx- 
imately 2 GeV. One way to compensate for this effect 
is to shift the px range at ^/s = 630 upwards by a few 
GeV. Due to the steep negative slope of the jet-pT spec- 
trum, it is sufficient to shift only the lower edge of the pT 
bins. When this is done, Fig. [l2| shows that the change in 
gluon-jet fraction is A/ < 0.03. We do not correct / for 
this, but account for this residual effect in the systematic 
uncertainty associated with the jet px- 

Changing the gluon-jet fractions used in the analysis 
gives a direct estimate of the uncertainty on the subjet 



FIG. 28. The subjet multiplicity in Monte Carlo quark 
jets. The extracted distribution has been unsmeared. The 
tagged distribution was obtained directly from particle-level 
quark jets at ^ = 1800 GeV. 
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FIG. 29. Corrected subjet multiplicity for gluon and quark 
jets, extracted from D0 data. 



multiplicity. We will motivate the range of uncertainty in 
gluon-jet fractions at the two center-of-mass energies by 
investigating the behavior of the PDFs. For the jet sam- 
ples used in this analysis, the average jet px was approx- 
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PDF set 


(GeV) 


X 


xg{x) 


e 




CTEQ4M 


1800 


0.07 


1.643 


0.00 


0.59 


CTEQ4HJ 


1800 


0.07 


1.643 


0.00 


0.59 


CTEQ2M 


1800 


0.07 


1.714 


0.04 


0.60 


CTEQ5M 


1800 


0.07 


1.614 


-0.02 


0.59 


CTEQ5HJ 


1800 


0.07 


1.586 


-0.04 


0.58 


MRST5 


1800 


0.07 


1.586 


-0.04 


0.58 


GRV94 


1800 


0.07 


1.743 


0.06 


0.60 


CTEQ4M 


630 


0.2 


0.365 


0.00 


0.33 


CTEQ4HJ 


630 


0.2 


0.340 


-0.06 


0.32 




630 


0.2 


0.385 


0.06 


0.34 


CTEQ5M 


630 


0.2 


0.340 


-0.06 


0.32 


CTEQ5HJ 


630 


0.2 


0.350 


-0.03 


0.32 


MRST5 


630 


0.2 


0.290 


-0.21 


0.28 


GRV94 


630 


0.2 


0.405 


0.12 


0.36 



TABLE II. Values of gluon-jet fractions for different 
PDFs, calculated using Eq. ( ^.2[ ), at a jet pr = 65 GeV. 
The CTEQ4M parameterization is chosen as the reference. 
The fractional change in the gluon PDF g{x) is given by 

where g^°^(x) is the reference. 



imately 65 GeV. This jet px probes an average x value 
of 0.07 at Vs = 1800 GeV and 0.2 at = 630 GeV. In 
these regions of x, the quark PDFs are well-constrained 
by existing data. However, the gluon PDF is not so well- 
constrained. We examined different parameterizations of 
the gluon PDF at the two x values of interest. In par- 
ticular, the MRST5 gluon PDF is 21% smaller than 
the CTEQ4M parameterization at x = 0.2, but only 4% 
smaller ai x — 0.07. This and other comparisons between 
PDFs show larger fractional differences aX x = 0.2 than 
at X = 0.07. 

Assuming that the quark distributions are essentially 
identical in different PDF parameterizations, the gluon- 
jet fraction / for different PDFs can be estimated as 



/ 



rcf 



(Jref + e/-f) + (l-/'-°0 



(5.2) 



where /''''^ is the gluon-jet fraction from some reference 
PDF, and e is a fractional difference in the gluon PDF. 
Table H shows the gluon-jet fractions estimated for PDFs 
at the two center-of-mass energies. The MRST5 set 
shows the largest departure relative to CTEQ4M. In all 
cases, the change in / is in the same direction at both 

The preceeding discussion assumed that the PDFs 
had the same quark distribution. In reality, the quark 
PDFs also tend to change when the gluon PDF changes. 
Whe n th is compensating effect is taken into account in 
Eq. (5_^), the equivalent MRST5 gluon-jet fractions be- 
come /i8oo = 0.58 and /eao = 0.29. 

Based on the above, we assign uncertainties to the 
gluon-jet fractions of ±0.02 at = 1800, and ±0.03 
at -^5 = 630. In fact, we vary the gluon-jet fraction in 
opposite directions, using /isoo = 0.61 and /eao = 0.30, 
and /i8oo = 0.57 and /eso = 0.36, to gauge the impact on 



/l800 


/eao 


{Mo) 




r 


0.59 


0.33 


2 21 -1- n 03 


1 69 + n 04 


1 75 + n 1 5 


0.61 


0.30 


2.18 ± 0.02 


1.72 ± 0.04 


1.65 ± 0.12 


0.61 


0.36 


2.20 ± 0.03 


1.67 ± 0.05 


1.79 ± 0.17 


0.57 


0.30 


2.21 ± 0.03 


1.70 ± 0.04 


1.72 ± 0.14 


0.57 


0.36 


2.24 ± 0.04 


1.65 ± 0.05 


1.92 ± 0.22 



TABLE III. Subjet multiplicity in gluon and quark jets, 
and their ratio, extracted from D0 data and corrected to the 
particle level, assuming different gluon-jet fractions at the two 
center-of-mass energies. 



r. As in Sec. V A , we repeat the analysis assuming these 
different input gluon-jet fractions, this time including the 
correction to the particle level. The extracted ratios are 

The largest departures from 
1.75 define the systematic un- 



summarized in Table III 
the reference value of r = 



certainties of ±o;io- 

The second-largest source of systematic uncertainty in 
the subjet multiplicity stems from an uncertainty in the 
measurement of jet pt- A mismeasurement of jet pr will 
lead to the selection of a slightly different sample of jets, 
but will not affect the subjet multiplicity directly. If jet 
Pt is mismeasured at both center-of-mass energies, we 
expect the effect to partially cancel in the ratio r. An 
estimate of the impact from this uncertainty is therefore 
obtained by varying the jet px only at ^/s = 1800 GeV. 
Since the calorimeter response is independent of -y/s, we 
estimate the effect of a difference in any offset in px at 
the two center-of-mass energies by changing the jet-pT 
window from 55 < pt < 100 GeV to 57 < < 100 
GeV at = 1800 GeV. A 2 GeV shift in the measured 
jet Pt corresponds approximately to two times the total 
offset po for k± jets reconstructed with D — 0.5. [This 
assumes po{D) scales as D^po{D = 1.0)]. This reduces 
the subjet multiplicity ratio r by 0.12, which is taken as 
a symmetric systematic uncertainty. 

Because the correction to the particle level produces a 
large change in the shape of the subjet multiplicity dis- 
tribution, we will estimate the impact of the unsmearing 
on the systematic uncertainty on r. This uncertainty 
has two parts: one is the uncertainty due to the simu- 
lation of effects arising from dependence on luminosity, 
and the other is the uncertainty in the simulation of the 
D0 calorimeter. To account for the former, we use an 
alternate Monte Carlo sample at y/s = 630 GeV, with 
a luminosity of £ « 0.1 x lO'^^cm^^s^^, and note that 
r increases by 0.13. Such a small change in r indicates 
that it depends only weakly on luminosity. Nevertheless, 
we increase our nominal value of r = 1.75 by half of the 
difference (to r — 1.82), and take this correction as a 
symmetric systematic uncertainty of ± 0.07. 

To evaluate the other part of the uncertainty on the 
unsmearing, we compare two types of simulations of the 
D0 calorimeter. The default fast simulation (SHOWER- 
LIB) is a library that contains single-particle calorimeter 
showers obtained using the GEANT full detector Simula- 
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Source 


5r 


Gluon-jet fraction 


+U.i7 
—0.10 


Cutoff on jet pr 


±0.12 


Unsmearing 


±0.07 


Detector simulation 


±0.02 


Total 


+i).'2-2 
-0.18 



TABLE IV. Systematic uncertainties on the ratio r. 



tion. SHOWERLIB truncates the number of calorimeter 
cells associated with each individual particle, but rescales 
the energy of the shower to agree with the average energy 
given by the full GEANT simulation. The full GEANT 
simulation, while slower, accounts for the precise geome- 
try of the uranium plates in the calorimeter and has no 
truncation. In a test using a limited number of Monte 
Carlo events, the latter simulation produced more subjets 
than the former, and so we increase the value of the ratio 
by 0.02 (half the difference of the r values in each simu- 
lation) to r = 1.84, and take this correction as another 
symmetric systematic error of ± 0.02. Applying the same 
additional corrections to the nominal ratio in the Monte 
Carlo gives a final result of r = 1.91 for HERWIG. 

A list of the systematic uncertainties is shown in Ta- 
ble rV , all of which are added in quadrature to obtain the 
total uncertainty of the corrected ratio. The final result 
for the ratio is 



^ (Mg) - 1 

(M,) - 1 



1.84 ±0.15 (stat.) ±^;f| (sys.). (5.3) 



VI. CONCLUSION 

We present two analyses of D0 data using the k± jet 
reconstruction algorithm. One analysis examines the px 
and direction of k± jets reconstructed with the parameter 
D = 1.0. For this measurement of the jet px spectrum, 
we describe a procedure to calibrate the momentum of kj_ 
jets based on our experience with the cone algorithm, but 
using an improved technique for determining the offset 
correction. Compared to our published results for the 
cone algorithm with TZ = 0.7 the k± jet algorithm 
with D = 1.0 reconstructs 40 — 50% more energy from 
uranium noise, pile-up, multiple pp interactions, and the 
underlying event, and has a smaller uncertainty on the 
offset. We also report the results of a direct comparison of 
the k± and cone algorithms, on an event-by-event basis. 
Considering only the two leading jets in the central region 
{\ri\ < 0.5), the k± and cone jet axes coincide within 
ATI = 0.1 (0.5) at the 91% (99.94%) level. Matching 
with ATZ = 0.5, the corrected pt of k± jets is higher 
than the corrected Et of cone jets. The difference is 
roughly linear in jet px, varying from about 5 GeV at 
PT « 90 GeV to about 8 GeV at pt « 240 GeV. 

In the other analysis, we probe the structure of central 
k± jets reconstructed with the parameter D = 0.5, and 



find that the HERWIG Monte Carlo predictions of sub- 
jet multiplicity are in excellent agreement with our mea- 
surements. The subjet multiplicities in gluon and quark 
jets, predicted by a fully resummed calculation fl^ , and 
shown in Fig. BOl are qualitatively consistent with our 
data, but their mean values are slightly high. This dis- 
crepancy may be due to the fact that the calculation 
lacks a preclustering algorithm. The subjet multiplicity 
distributions, where we have subtracted the D0 values 
from the predictions, are shown in Fig. ^ The ratio of 
mean multiplicities for the resummed calculation (which 
assumes M < 5) is r = 2.12. The ratio in the D0 data 
increases by 0.06 with the assumption M < 5. There- 
fore, the resummed prediction is well within the limits 
of experimental uncertainty. The ratio measured at D0 
agrees with the result of r = 1.7 ± 0.1 from ALEPH, 
measured in e'^e~ annihilations at ^/s ~ Mz for a subjet 
resolution parameter yo — 10""^ ||l^, and with the associ- 
ated Monte Carlo and resummation prediction |14| , but 
is higher than the ratio measured at DELPHI ]1(]| ]. The 
DELPHI result uses a different definition of the jet reso- 
lution scale y than used by ALEPH (yi ) , which takes the 
place of Z? in a hadron collider, making direct compar- 
isons difficult. These experimental and theoretical values 
for r are all smaller than the naive QCD prediction of 
the ratio of color charges of 2.25. This may be caused 
by higher-order radiation in QCD, which tends to reduce 
the ratio from the naive value. 

In summary, we present the first detailed measure- 
ments of properties of /cj^ jets in hadron collisions. Using 
the standard value D — 1.0 of the jet-separation param- 
eter in the k^ algorithm, we find that the pt of fcj^ jets is 
higher than the Et of matched cone jets (with Tl, — 0.7) 
by about 5 (8) GeV at « 90 (240) GeV. To analyze 
internal jet structure, we measure the multiplicity distri- 
bution of subjets in k\_ jets with D = 0.5 at ^/i = 1800 
GeV and 630 GeV. Exploiting the difference in gluon-jet 
fractions at the two center-of-mass energies, we extract 
the subjet multiplicity in gluon and quark jets. The dif- 
ferences between gluon and quark jets are summarized in 
the ratio of average emitted subjet multiplicities, mea- 
sured as: 



1 



1 



1.84 ±: 



0.27 
0.23 



(6.1) 



The D0 result demonstrates that gluon and quark jets 
are significantly different in hadron collisions, and that 
it may be possible to discriminate between them on an 
individual basis. 
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FIG. 30. The subjet multiplicity in g luon and quark jets, 
for i/cut = 10~^ (as defined by Eq. [p.q]) , in a resummation 
calculation by Forshaw and Seymour ]I7|. The jets are pro- 
duced at = 1800 GeV, with pr = 65 GeV and = 0, using 
the CTEQ4M PDF, and are reconstructed with D = 0.5. The 
points in the fifth bin refer to M > 5. 
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FIG. 31. The subjet multiplicity in (a) gluon and (b) 
quark jets, for D0 data, for the HERWIG Monte Carlo, and 
resummed predictions. The resummed prediction does not 
use a preclustering algorithm. The points in the fifth bin are 
for M > 5. The D0 data (see Fig. |2^) have been subtracted 
from each set of points. 
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